The sinusoid and the phasor 
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Mathieu equation is widely used to study several natural phenomenon. In this paper, we show 
that replacing the sinusoid in the Mathieu equation with a phasor can lead to solutions that behave 
in a totally different way. Solutions of Mathieu equation are either bounded or grow unboundedly 
at an exponential rate. Solutions of this new equation are always unbounded and grow linearly with 
time. 

Linear differential equations with periodic coefficients have applications in many areas like Plasma 
Physics P, 0| and Condensed Matter Physics [3]. The general form for such an equation is 



X = A(t)x 



(1) 



where x is a vector and A(t) is an n x matrix with time-periodic coefficients. For n — 2, a, very 
important example of such a system of equations is the Mathieu equation 0, 



X — 2qxcosujt 



(2) 



which can be written in the matrix form as 

d_ / xi 

dt \ X2 



1 

2q cos ijjt 



Xi 
X2 



(3) 



where q is an arbitrary constant and uj is the oscillation frequency of the field. 

In many problems of interest in physics and engineering, the sinusoid is replaced by a phasor since it 
often makes the analysis easier. If we replace the sinusoid in Eq. © by a phasor, we get 



(4) 



where y is a complex function of a real argument, t. It must be noted that Eq. Q can also be written 
in form of Eq. ^ with n = 4. Transforming to the variable, t — J ^e^™*/^, Eq. Q becomes 



(Py dy 



(5) 



which is the well known Bessel equation The transformation, r = y^e-'"'*^^, used to get Eq. ^ is 
not one-to-one, but as we will see later, the solutions have a branch-cut, which permits the use of this 
transformation leading to non-periodic solutions. Equation ^ has two linearly independent solutions 
and any general solution can be written as a linear combination of these. 



y (r) = acj) (r) -|- 6^ [t) 



(6) 



where 0, -0 are Bessel functions. 



0(t) = J, 




]wt/2 



and V (''') = 




]wt/2 



(7) 



The velocity of the particle at any instant of time is given by, 

v{t) = a0'(t) + biP'{t) 



(8) 
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where, ' represents differentiation w.r.t. t . 

In order to find the path of the particles in phase space, we need to evaluate the bessel functions 
and their derivatives. The derivative of Jo is — Ji and that of Yo is —Yi. Though all these functions 
are analytic for complex arguments, Yq and Yi have a branch cut. Yq,Yi can be evaluated using the 
expression Q 



Yoir) 



(c + log0Jo(r)+2f;(-ir-i:^ 



n=l 



(9) 



and 



Yi (r) = 



C + log ^) Ji (r) - - Jo (r) + ('^2„-i (r) - J2„+i (r)) 



(10) 



where, C = 0.577215665 is the Euler's constant . Although Yq (r) and Yi (r) seem to be a periodic func- 
tions of t, the condition that y and v must be continuous and differentiable imposes certain restrictions. 
We need to take special care of the branch cut in log ^ term. In our problem, 



which implies 




(11) 



T Uq jwt 

log- = logW- + — 



(12) 



From Eq. (fT^ . it can be seen that the continuity and differentiability condition forces ■0(i) to be growing 
in time. Equation (|12p predicts this unbounded growth with time to be linear, whereas unbounded 
solutions of the Mathieu equation, Eq. ©, would have exponential-in-time growth. Also, Eq. Q has 
no bounded solutions for any value of q, whereas all solutions of Eq. ^ are bounded if \q\ is below a 
certain critical threshold. 

One of the authors (KS) would like to acknowledge support of the Israel Science Foundation and the 
Minerva Foundation. 



[1] W. Paul, Rev. Mod. Pliys. 62, 531 (1990) 

[2] K. Sliah and H. S. Ramachandran, Phys. Plasmas 15, 062303 (2008) 

[3] Charles Kittel, Introduction to Solid State Physics (Wiley, New York, 1996) 

[4] A. Abramowitz and A. Stegun, Handbook of Mathematical Functions (Dover, New York, 1964) 

[5] N.W. McLachlan, Theory and Application of Mathieu Functions (Dover, New York, 1964) 

[6] M. C. Gray, Comm. ACM 4, 169 (1961) 



